Fano interference and cross-section fluctuations in molecular photodissociation 
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We derive an expression for the total photodissociation cross section of a molecule incorporating 
both indirect processes that proceed through excited resonances, and direct processes. We show that 
this cross section exhibits generalized Beutler-Fano line shapes in the limit of isolated resonances. 
Assuming that the closed system can be modeled by random matrix theory, we derive the statistical 
properties of the photodissociation cross section and find that they are significantly affected by 
the direct processes. We identify a unique signature of the direct processes in the cross-section 
distribution in the limit of isolated resonances. 



Spectral correlations of closed quantum systems, whose 
associated classical dynamics are chaotic, are known to 
be nearly universal and can be modeled by the Gaussian 
invariant ensembles of random matrix theory 0, 0, 
When such systems become open through their coupling 
to continuum channels, their bound states acquire de- 
cay widths and become resonances, but they are still ex- 
pected to exhibit universal statistics 4]. Examples are 
the conductance fluctuations in quantum dots and 
the statistics of the indirect molecular photodissociation 
cross section 0,0] ■ If the coupling is weak, the corre- 
sponding resonances are isolated and are often charac- 
terized by a Lorentzian line shape (in quantum dots it is 
necessary to assume temperatures that are much smaller 
than the decay width). The statistics of the resonance 
widths are determined by the corresponding statistics of 
the bound states and energies £'„ of the closed system. 
For a classically chaotic system, these statistics can be de- 
rived from random matrix theory. In recent years, such 
a random matrix approach has been successfully used 
to model the statistics of resonances -^l and cross-section 
fluctuations in the pho todissociation of classically chaotic 
molecules [^l9[ lin[lll[| l ^ ll3j. A semiclassical treatment 
was discussed in R.efs.ll4l andll5[ 

However, as first observed by Beutler 0| and later in- 
terpreted by Fano 0, the line shape of an individual 
resonance may differ substantially from a Lorentzian: in- 
terference between the indirect decay via the quasi-bound 
state and direct (fast) decay to the continuum gives rise 
to a so-called Beutler-Fano line shape. For real bound- 
state wave functions, the line shape (versus energy E) is 
proportional to (g„ + e„)^/(l 4- e\) where e„ = 2{E — 
En)/^n- Here, g„ is the Fano parameter characterizing 
the line shape, and r„ is the width of the n-th resonance. 
Beutler-Fano profiles have been observed in molecular 
photodissociation [l^l, autoionization conductance 
through quantum wells [l^l and quantum dots [iJl , STM 
spectroscopy of surface states , semiconductor super- 
lattices [23 , and Aharonov-Bohm rings . 

The Fano parameter g„ fluctuates from one resonance 
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FIG. I: Indirect molecular photodissociation. Shown are 
three electronic surfaces: the ground-state surface eo (in 
which the ground state \g) resides), a binding surface ei and 
a repulsive surface 62. The surfaces ei and 62 give rise to 
an effective electronic surface with a barrier near their cross- 
ing (not shown), and the bound vibrational states |n) in 
ei become resonances. Indirect photodissociation proceeds 
through these resonances with a total cross section cy{E) 
(shown on the left) that is a sum of Lorentzian line shapes. 



to another. Distributions of the Fano parameter have 
been calculated for transmission through a quantum 
dot and for the photodissociation of molecules [l^ 
in the one-channel case (assuming that the corresponding 
closed system is classically chaotic). 

Here we derive an expression for the total photodisso- 
ciation cross section for any number of open channels in 
the presence of direct decay processes (see Figs. Q] and |21l. 
We show that in the limit of isolated resonances their line 
shapes have the form of generalized Beutler-Fano pro files 
+ £nP/(l + £^) with a compZex Fano parameter |25j qn- 

Given that direct photodissociation processes affect 
the line shapes so dramatically, it is interesting to find 
out how they affect the statistics of the photodissocia- 
tion cross-section when the closed system is classically 
chaotic. We derive a closed expression for the cross- 
section autocorrelation function (in energy) and find 
that it is universal provided that the excitation process 
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and the continuum coupling are spatially well-separated. 
System-specific information enters only in the values of 
the direct and indirect channel couplings. We also calcu- 
late the cross-section distribution, and find that the di- 
rect decay gives rise to a characteristic maximum in the 
distribution (see Fig. ^ in the regime of isolated reso- 
nances. This is in contrast to the monotonically decreas- 
ing behavior of the distribution in the same regime, in 
the absence of direct coupling. In summary, we show that 
cross-section fluctuations are significantly affected by the 
presence of direct decay channels, but remain universal. 
Our results also apply to atomic autoionization |0, . 

A molecule can dissociate into several channels c by 
absorbing a photon. In the dipole approximation, the 
total photodissociation cross section at energy E is given 
by 



a{E) = ME)J2\{<^i-\Em9)\' 



(1) 



where fl — p,-e is the component of the dipole moment fi 
of the molecule along the polarization e of the absorbed 
light and (Tq{E) cx {E — Eg). Here \g) is the ground state 

with energy Eg, and |$c \e)) (c = 1, . . . , A) is a disso- 
ciation solution at energy E defined by an outgoing wave 
in channel c and incoming waves in all other channels. 
We consider a model 0, la, S 113 in which the Hilbert 
space is divided into two parts: an internal "interacting" 
region, and an external "channel" region (cf. Fig.|2Il. The 
internal region is described by the Hamiltonian Hq rep- 
resented by an TV X matrix Hq with cigcnstates \n) 
(n = 1, . . . , A^). The external region is spanned by the 
A open dissociation channels |c). The two regions are 
coupled by an operator W that can be represented by an 
A^x A matrix W with matrix elements (c|Vl^|n) = jnc- In 
general, the dipole operator p, can couple the ground state 
to both the internal states \n) and the external channels 
|c). We define \a) = fJ,\g), and introduce two vectors a'" 
and a"^^. The first has A^ components a™ = {n\a), de- 
scribing the dipole coupling to the internal states, and 
the second has A components af^ = {c\a), describing the 
dipole coupling to the continuum channels. 

We first show that, in the regime of isolated reso- 
nances and for a'^'^ ^ 0, the cross section describes 
a sum over Beutler-Fano resonances. We write 
'E^^,\<i>i~\E)){<^i-\E)\ = -Tr-nmG{E + iO) and sep- 
arate the channel and internal components of the Green 
hmction G{E + iO) We obtain 



aiE)/aoiE) Im[(a|Gch|«) 

TT 
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FIG. 2; Top: states and transition matrix elements in the 
random matrix model discussed in the text. Bottom: the 
photodissociation cross section (j{E) (arb. un.) for one real- 
ization of Hq (circles) and for (5 = 1, N = 128 and A = 1. 
The solid lines describe fitted Beutler-Fano liner shapes. 



structured decay continua, Eq. ||2Jl simplifies to 

a{E)/ao{E) = 



(3) 
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Here H^s ~ Hq — ittWW^ is an effective (non- 
Hermitean) N x N Hamiltonian in the internal space. 
In the absence of direct photodissociation, a."^^ = 0, and 
Eq. ||2Jl reduces to the result of Refs. and0. 

In general, the resolvent in Eq. Q can be written 
as {E - H,s)-^ = J2n \Rn){Ln\iE - £n)-\ whcrc |i?„) 
and |L„) arc bi-orthonormal right and left eigenvectors 
of ffoff with complex eigenvalues £n- However, in the 
regime of isolated resonances we can apply the Breit- 
Wigner approximation |i?„) « |L„) « |n) and Imf„ = 
— r„/2 « — X)c l7ncP- The cross-section Q can then 
be written as a sum over resonances. In the presence 
of direct photodissociation, the contribution from each 
of these resonances has a generalized Beutler-Fano line 
shape |g„ -I- e„p/(l -I- e^) with a complex parameter g„ 
whose real part and modulus are given by 



Re q„ 



„ch|2 



(4) 



where Gch is the channel Green function. Assuming un- 



In general, Im qn ^ 0, and there is no energy for which 
the cross section vanishes. However, for p — I and 
A = 1, Eq. Q simplifies to Fano's expression Req„ — 
a™/{naf^jnc) and Imqn — (assuming that all matrix 
elements are real). 

In the following we calculate the statistical proper- 
ties of the photodissociation cross section assum- 
ing that the dynamics in the closed interaction re- 
gion are fully chaotic. The matrix Hq is taken to be 
a. N X N random matrix from the Gaussian Or- 
thogonal Ensemble (GOE) or from the Gaussian Uni- 
tary Ensemble (GUE), with distribution P{Hq) dHo cx 
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exp[-(/3iV/4)Trff2]dHo- Here /3 = 1 in the GOE and 
P = 2 in the GUE. In the hmit of large N, the aver- 
age eigenvalue distribution (normalized to 1) is i^{E) — 
(27r)~^V4 — E'^ for \E\ < 2 (and zero otherwise), and 
the corresponding mean eigenvalue spacing is A{E) ~ 
l/[Ni'{E)]. In the same limit (i.e., for large TV), a™ and 
7nc (c — 1 , . . . , A) are independently distributed Gaus- 
sian random variables, and for each n ISl: 



^«">7n)oc exp 



/3 



(5) 



Here7t = (^ti, • • ■ , T^a) and 

M = N-^ V^V with y = (a'", W) (6) 

is an (A 1) X (A + 1) matrix. In the following we as- 
sume the channel vectors (columns of W) to be mutually 
orthogonal. This can always be achieved by a suitable or- 
thogonal (unitary) transformation in channel space. 

Average cross section. In the center of the band {E = 
0), the average cross section is (ct) = co [ I P/'''' + 

Ec=i |a?lV(l + Ac)] = <Jind + Ec'^c", where aind and 
a'^" are the average cross sections in the limiting cases 
of purely indirect and purely direct dissociation, respec- 
tively. Here Ai, . . . , Aa are the A dimensionless eigenval- 
ues of the matrix ■e^v'W'^'W . It is often convenient to 
characterize the strength of the coupling to the contin- 
uum by transmission coefficients = 4Ac/ (1 -I- Ac)^. 

Cross-section autocorrelation function. We define _a di- 
mensionless cross section autocorrelation function 



S{E,uj)^ar^l[{a{E-u;/2)a{E+uj/2))-{a)^] . (7) 

In the Breit-Wigner approximation, this correlation is 
most conveniently calculated in the time domain js^]. 
Defining C{E,t) = /^dwe''^* S{E,uj), and using © we 
find 
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C{E,t) = — [Apit)-Bj{t)b,^p{t)] 



(8) 
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FIG. 3: Cross-section autocorrelation functions. The sym- 
bols are from simulations of random matrices {N = 128) 
at the center of the band {E — 0) and for A = 10 and 
Ac = 10"'^ for all c. The autocorrelation function without 
a direct coupling (open squares, 8 — 0), and with direct cou- 
pling {6 — 1.26, open circles) are compared to Eq. 11211 (solid 
lines) . 



where &2,/3(i) is the two-level form factor 0, and Ap{t), 
Bp (t) are functions that depend in general on the matrix 
M in Eq. Q and the dipolc coupling coefficients to the 
continuum af^. 

The expression for Ap{t) and B0{t) simplify when the 
dipole "channel" is orthogonal to all channel vectors, i.e., 
TytQ,in ^ This is the case when the excitation process 
and the continuum coupling are spatially well separated. 
Using the so-called rescaled Breit-Wigner approximation 
0,113, we find 



Mt) - Y[{l + 2T,t)-'/^h + -J2rcil + 2T,t)-' 



^(5:r,(l + 2T,i)- 



Bi{t) = l[{l + T,t)-'/'[l--Y,rc{l+T,t)-^ 

c—l c 

for /3 = 1 . A similar result is obtained for /3 = 2 

A 

Mt) = Y[{l+T,t)-'[2+-(^Y.^,{l+T,t)-' 



(9) 



c=l 
A 



(10) 



Mt) = Y[i^+Tct/2)-'[l--Y,^c{l + T,t/2)-^ 

c—l c 

The results H9I10|I describe universal correlations; they 
depend on the coefficients Tc and on the parameters 
Tc/Tc = crf^/cTind which measure the strength of the di- 
rect photodissociation channels, but they do not depend 
on the microscopic details of the system (such as the 
ground state or the nature of the excitation mechanism). 

Two special cases of (O and HlOfl are of interest. First, 
in the limit of Tc = 0, the results of Ref. are recovered, 
valid in the absence of direct coupling to the continuum. 
Second, consider the case of A equivalent open channels 
= T, erf = a'^'\ In the limit A oo , T ^ with 
AT = K constant, the expressions in Eq. jnj simplify to 



Ap{t) = Ape-''^'^ , Bp{t)^Be 

with = 3 + 6i/2 + 36'2, A2 = 2(1-^0/16), B 
and 9 ^ At = na'^" /aind- We obtain 



|t|/2 



(11) 



(1-0/4), 



S{E,io) = ^ Apfiuj) 



(12) 



B' 



d( / 
TT 



where Y2,f3{uj) is the two- level cluster function T|, and 
ficu) = {k/2)/{uj^ + K^/A). Fig.Hshows Eq. ^ (solid 
lines) together with results from random matrix simula- 
tions (symbols). For /? = 1 and in the presence of direct 
decay, the correlation function is close to a Lorentzian. 

Cross-section distribution. The distribution P(^a/aind) 
is calculated from its Fourier transform -F)3(s) = 
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^g-isa/ai„d^ within the Breit-Wi gner approximation. We 
have calculated ^/3(s) for A equivalent open channels in 
the limit of A — > cx) with AT = k kept constant. In this 
case, r„ ~ r = 2AX/N, and using Q we obtain 



^gise/(2rw) 



det[{E-Hoy 




(13) 



with r| 



r + 47ris/(iV/3))(r-27ris6'/(jV^)). Eq. d 
can be evaluated using the results of Ref. |33- For [3 — 2 



F2{s) 



gise/(2rAf)g-7rr2/A 



cosh(f)+isinh(^)(i- + 



r 



12 

r 



.(14) 



For (3 = 1 the corresponding result can be expressed in 
terms of a four- fold integral [s^- Fig. 0] shows the in- 
verse Fourier transform of Ffj{s) (solid lines) in compar- 
ision with random matrix simulations (symbols). In the 
presence of direct coupling, the cross-section distribution 
exhibits a maximum (see Fig. ^Jj). In the limit of iso- 
lated resonances, this maximum is a clear signature of 
the direct processes. 

In conclusion, we have shown that a direct coupling to 
the continuum leads to generalized Fano resonances in 
the total photodissociation cross section, and used ran- 
dom matrix theory to derive the signatures of these direct 
processes in the cross-section statistics. 
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FIG. 4: Cross section distributions, (a) The cross section vs. 
energy E in the absence of direct coupling {8 — 0) and for one 
random matrix realization of Ho- Here (3 = 1, N = 64, A — 
10, and Ac = 5 X IQ-^^ for all c. (b) Same as in (a) but in 
the presence of direct processes {6 = 0.125). (c) Cross-section 
distributions at i5 = for 9 = (squares) and 6 = 0.125 
(circles), (d) as in (c) but for (3 = 2. The solid lines in (c) 
and (d) are the inverse Fourier transform of Ff3{s). 
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